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SUMMARY 

The method presented in NACA TN 2692 for evaluating the skin fric- 
tion of a turbulent boundary layer in compressible flow on an insulated 
surface is extended to evaluate the turbulent skin friction and heat 
ransfer in compressible flow on a surface which is heated or cooled. 

The results of this analysis are in good agreement with the heat 

transfers measured in flight on the NACA RM-10 missile up to Mach num- 
bers 01 J.f. 


INTRODUCTION 


Within the last decade the tremendous increase in flight Mach num- 
bers that has been achieved has placed the problem of aerodynamic heating 
among the most important of the present time. In particular the problem 
o he heat transfer through turbulent boundary layers at high speeds has 
received considerable attention. It has long been known that the heat 
transfer through turbulent boundary layers is so intimately connected with 

he skin-friction force exerted by the boundary layer that the two problems 
are almost one. 


_. T of the turtulent boundary layer presented in NACA 

IN 2092 l ref. 1) permits an easy extension of the turbulent skin-friction 
law from incompressible to compressible flow. The analysis led to expres- 
sions for the extent of the effective laminar sublayer 5 ^ and for the 

velocity u L at the point S L . These results permitted the skin-friction 

law to be derived. If use is made of the quadratic dependence of temper- 
ature on velocity derived by Crocco in reference 2, the results of ref- 
erence 1 may be extended to obtain the heat transfer through the turbulent 
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boun dar y layer in both incompressible and compressible flows. It is 
the purpose of the present short analysis to make this extension. 


SYMBOLS 


parameter. 


n 


n(r - l) Vol 


1 

n+1 


.2 u-,5 


o 8 J 


c f 


skin-friction coefficient. 


2t 


w 


Po u o 



K 

k,m,n,r 

M 

N x 


N 6 

q 


r .f . 




T 


constant (O.OV?) 
constants 
Mach number 

qx 

Nusselt’s number, — — 

kAT 

c? 

Nusselt’s number, 

k AT 

local rate of heat transfer 
recovery factor 

Reynolds number, 

Reynolds number, 
absolute temperature 


u 0 x 

v o 

UqS 



NACA RM L52H04 


3 


u velocity in x-direction 

x distance along surface 

y distance normal to surface 

7 ratio of specific heats 

6 boundary- layer thickness 

6 boundary- layer momentum thickness 

k thermal conductivity 

[i viscosity 

V kinematic viscosity 

p density 

r shear stress 

Subscripts: 

av average value 

adw adiabatic wall conditions 

L , conditions at edge of laminar sublayer 

o free-stream conditions 

w wall conditions 

DERIVATION OF THE HEAT -TRANSFER LAW FOR TURBULENT BOUNDARY LAYERS 

Before the present analysis is presented, a summary of the findings 
of reference 1 is in order • It was found that for boundary layers of 
the type 
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the dimensionless effective extent of the laminar sublayer 
given by 



was 


£l 

6 


n(r - 1) V L 
k 2 u Q S 


n 

n+1 


and the velocity at the point was given hy 


hi = fn(r - 1) v 'l1 

u o L k2 u o 5 _ 


1 

n+1 


( 2 ) 


( 3 ) 


With these results, it was reasoned that the skin friction was 


w 






which, when the necessary substitutions had been made, resulted in the 
following skin-friction formula 


c f 


2 




2 

n+i 

Pq 



2 

n+1 


( 5 ) 


Since 


PL = ^o 
p o 


(6) 


and if it is assumed that 



( 7 ) 
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then 


Cf = 2 


n(r - 1) 


1-n 

n+1 


% 


2 n-gm-1 

n+1 /T 0 \ n+1 

W 


( 8 ) 


For a velocity profile having n — 7* it was found in reference 1 (from 
skin-friction measurements in incompressible flows) that 


n(r - 1) 
k 2 


158 


so that the local skin-friction coefficient for the compressible case 
became 


-1/k 

Cf = 0.0^5 Rg 



( 9 ) 


With these results it is now possible to pass on to the derivation of 
the heat transfer through the turbulent boundary layer. 

If a temperature distribution in the boundary layer similar to 
that used by Crocco (ref. 2) for Prandtl number = 1 and zero pressure 
gradient is adopted; namely: 


the temperature in the boundary layer may then he written 


T = T w - (t w - T adw )^ - (T adw - T 0 ) ^) 2 (10) 

It might be noted that, although the quadratic form of the tempera- 
ture dependence on velocity was derived under conditions of zero pressure 
gradient, this form appears also to be a very useful relationship for 
most supersonic-missile shapes where the pressure gradients are generally 
quite small. J 
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2 - = 0 


Differentiating equation (10) and evaluating it at the wall where 
gives 


St\ _ _ / T w ^adw\ /du\ 
& /w \ u 0 /Ww 


( 11 ) 


Hence, the heat transfer at the wall is 


~%i ~ K w 



K w fry ~ I'adw) /Su\ 

u 0 \Sy/w 


( 12 ) 


From equation (4) it will be seen that 

(b u\ = ^L U L 
\6y/ w ^ W&L 


(13) 


so that 


% 


_ k w^L 


^ 5 l 



m 


If the values of S T and — that are obtained: from equations (2) 

u o 

and (3) are substituted into equation (l4) and / the heat transfer is 
expressed in terms of a Nusselt number based on the dimension &, the 
result is 


N 6 


q& 

k cA t 


n(r - l) 
k 2 


1-n 

n+1 



n-1 

n+1 


i^w 

K o 


(15) 


If it is assumed that 


K w ^ Vu 
*o ~ ^o' 


(16) 
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equation ( 15 ) may be written 


N 5 = 


n(r - 1) 


IzIL n-1 
.1 n+1 


k2 




n+1 


n-1 

^\ n+1 ^ 

W Hn 


(17) 


n-1 
/ v \ n+1 

Substituting the value of ( v j 
equations (6) and (7) results in 


HL 


and — obtained from the use of 


1-n 


n-1 


N 5 = 


n(r - 1 

iT + 1 b 5 ct M 

k2 

J E5 t l 


n-2m-l 
n+1 


( 18 ) 


Writing the equations for dimensionless skin friction and heat 


transfer together and expressing 2 


n(r - 1) 


k2 


i 


1-n 

n+1 


as a constant yields 


_ 2 

CjF = Constant x R§ n+ ^ 

n-1 

jjgp _ Constant n+1 


(19) 

( 20 ) 


where 


n-2m-l 



Since the right-hand sides of equations (19) and (20) are the 
incompressible expressions for cp and N & , respectively, the effect 

of Mach number and temperature potential must be contained in the 
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factor F. This factor may be evaluated in the following manner. From 
equation (10) 


Tl T w (t w ^adw) U L _ (^adw ~ T Q ) / u l \^ 


T 0 T 0 T 0 

u o 

T 

± o 

which may be written 



— = 1 + r ' (7 - l)M^ 

1 - (if 

T 

+ w 

1-3 

0 

ro 

L v u °/ _ 



\u. 


- T 


adw I 2 . _ 

u,. 


(22) 


(23) 


The factor F therefore becomes 


F = < 1 + Z-JLl ( 7 - l)M S 


■ - sf 


n- 2m-l 
n+1 


' 5 ) f (24) 


It may be seen that the second term within the brace represents 
the effect of Mach number, and the third term the effect of temperature 
potential on the skin-friction and heat-transfer coefficients. The 

value of — is usually of the magnitude 0.4 and although it depends 


u r 


fp “ ^adw 

upon M, Rg, and — — , it generally does not have a first-order 


u r 


Ut 

effect upon the factor F. The value of however, must be found 

in the following way: 

Since 


Ur 


U r 


“1 


1 

n+1 


n(r - 1) VL 1 


u o 5 J 


(3) 
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1 1 
n+1 \n+l 


_ / nfr - i) jA m - A n+1 

u o \ k 2 R&j \ v o/ \ T o/ 


m+1 

\n+l 


(25) 


/ v L\ n+1 

where the value of 1—1 has been expressed in terms of temperatures 
by means of equations (6) and (7) and where 


k 2 R & 


(26) 


Tt 

Now, if the value of — given by equation (23) is substituted in 

To 

equation (25) there results 


m+1 

n+1 


1 

r 


'I 

^ = A^ < 

1 + r,f — (7 ~ T). m 2 

1 - ^ 

... T w - T ad „ / . “L\ 

uo 

2 

\ u o/ 

T ° l V 




J 


(27) 


Ut 

This equation for — may be expressed as 


n+1 

(M m+1 + A m+1 r.f . (y - 1) m 2/^l\ 2 + A m+1 T w ~ T adv u l 
lun/ 2 \ u o/ T 0 Uq 


1 

a m+1 li + 

2 


_l)_ m 2 + T w ~ T a d.w\ 
T 0 


= 0 


(28) 


U L 


and solved graphically for — 


Tw " T a( jw 


for given values of Rg, M, and 
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For the particular case of a turbulent boundary layer with a one- 
seventh power velocity profile in air when 7 = 1.4 and m may be 
taken as 0.7 6 (an approximation which is usually adequate for calcula- 
tions such as these) , the equations necessary for the calculation of 
heat transfer become 


0.75 , % 

N&F = 0.0225R & (29) 

where 



Ut 

and — is found from the following equation 
u o 


M 1 *- 55 

VU rs) 


* (if 568 * 


XI £/Vf + /l58\°- 568 Ty - Tadw /M _ 

5 \ u o/ / T 0 \u 0 / " 


^158 \°- 568 / + r.f.M 2 + T w - T afa \ = Q 

^6 j V 5 . T 0 / 


(31) 


Ut 


Tables I and II give values of — and 


u^ 


found by using 


equation (31) for Mach numbers up to 5 for a range of values of Rg 
Ttt - T fl 

and that are useful in heat- transfer calculations. 


COMPARISON WITH EXPERIMENT 


In order to compare the results of this analysis with experiment, 
it is necessary to have measurements of local rates of heat transfer 
for conditions of turbulent flow when the local values of R 5 , M, and 
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Tw " T a dw 

To 


are known. 


This type of information for a range of Mach 


numbers, Reynolds number, and temperature potentials can be obtained 
from the heat-transfer and skin-friction measurements made on the NACA 
RM-10 missile in flight and reported in part in references 3 and 4. For 
two of the body .stations at which local heat transfers were measured in 
reference 3 (stations 85 and 122 inches from the missile nose) the bound- 
ary layer had been surveyed under similar conditions for the skin-friction 
study reported in reference 4, so that the boundary- layer thickness 6 
and hence R5 were known. The pertinent measured quantities at these 
stations for several flight conditions are given in table III . The first 
six points are taken from data published in references 3 and 4. The last 
four points are computed from heat- transfer data not yet published. The 
correlation of these data by the present method is shown in figure 1 
wherein the measured and theoretical results are plotted in the form N$F 
versus Rs. A reasonably good correlation of the data is obtained for 

T — T h 

Mach numbers ranging from 1.6 to 3-7 and temperature potentials — — 

T o 

ranging from 0.15 to -1.8. In comparing these experiments with the 
theory, it was assumed that the boundary- layer profile had a one-seventh 
power profile even though the surveys show that the power of the boundary- 
layer profile varied somewhat from test to test. This is not considered 
to be a serious matter in making a comparison between the theory and 
experiment, as experience has shown that small variations of profile power 
from the value of seven do not materially affect the magnitude of the heat 
transfers involved. 


In general, heat- transfer data are not plotted as has been done in 
figure 1, in terms of local correlations, but in terms of N x and R x . 

This usual practice is generally permissible in the case of a flat plate 
having no pressure gradient, but it may be of interest to see how the 
results of the present analysis appear when integrated along such a flat 
plate so as to be presented in more conventional form. The integrations 
necessary (carried out in detail in the appendix) result in the following 
relations for incompressible flow with n = 7 


c f = 0.0578 R x 


■1/5 


V 5 

N x = 0.0289R X 


Thus the normal minus one-fifth and four-fifth power variations of 
Cf and N x with R x are found. 
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For the case of compressible flow, the solution is not quite so 

1 ^ IJI n 

simple and only in the case of a flat plate where — — is con- 

T 0 

stant can the following useful approximations be made. For n = 7> 


-l/5 

c f G = 0.0906R X 

, 4/5 

N X G = 0.0453R x 


where 



It is, however, desirable for the sake of accuracy and generality 
to retain the relations given in equations ( 29 ), (30), and ( 31 ) and 
correlate turbulent-boundary-layer heat transfer and skin friction on 
the basis of the length . & rather than x. 


CONCLUSIONS 


The method presented in NACA TN 2692 for evaluating the skin friction 
of a turbulent boundary layer in compressible flow on an insulated surface 
is extended to evaluate the turbulent skin friction and heat transfer in 
compressible flow on a surface which is heated or cooled. 

The results of this analysis are in good agreement with the heat 
transfers measured in flight on the NACA RM-10 missile up to Mach numbers 

of 3 . 8 . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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• APPENDIX 

DERIVATION OF DEPENDENCE OF SKIN FRICTION AND 
HEAT TRANSFER ON R x FOR A FLAT PLATE 


The momentum equation for the boundary layer on a flat plate may 
he written 


d £ = Cf 
dx ~ 2 


(Al) 


Since 


£ 

6 


is a constant on a flat plate at a given Mach number 


d& = _5_ r 
dx 26 f 


Substituting the expression 


(A2) 


Cf 




into (A2) there results 


dS 

dx 


8_ k/v 0 

26 F \u 0 S 


2 

n+1 


(A3) 


(Al) 


When M = 0 and — = 0 , F 

T 0 


1 so that equation (Al) may be 


integrated to obtain 
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S 


n+1 n+1 




n+1 

n+3 


(A5) 


The local skin friction and Nusselt number may then be expressed for 
M = 0 and T w = T a dw as 


n+1 2 2 


- ■ : 5 ?)"(£)“ 

(A 6) 

n+1 -2_ n+1 

„ K n+3 /n + 1 2e N \ n+3 /uo*\ n+3 

■ -T\n ♦ 3 J (v 0 ) 

(A7) 

0.045, and ~ = -2- so that 
* & 72 


_ -l/5 
c f = 0.0578R X 

(A8) 

4/5 

N x = 0.0289R X 

(A9) 


These formulas are the more conventional expressions for local skin- 
friction coefficient and Nusselt numbers as a function of Reynolds 
number . 

It may be seen, from the differential equation (A3) that the extremely 
simple expressions just derived are not valid for the case when M or 

I |~^ ^ f | 1 -I 

^ is other than zero. However, F is not a very sensitive 

function of R$ and, if an average value of R& for a given problem 
is used to evaluate an F = F av , the resulting approximate equations 

are extremely useful and fairly accurate. Thus, for a flat plate with 
a constant surface temperature 
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and 


5 = 


n+1 



n + 3 

n + 1 



n+1 

n+3 



n+1 2 n+1 

= 1/ K \ n+3 /n + 1 2e\ n+3 /up*\ n+3 

2 Vav/ \n + 3 &/ \v 0 J 


(A10) 


(All) 


(A12) 


For n = 7 these expressions become 


c f (1)' 1/5 ( F av) = 0.0906 R x 1/5 

(A13) 

Nx (l) l/ 5 ( F av ) V5 = °- 0453R x 4/5 

(Al4) 

These equations may be written 


c f G = 0.0906 r x -1 / 5 

(A15) 

N X G = 0.0453R x 4 /5 

(A16) 


where 



(A17) 


In equation (A17) 


av 


U L 

is an average value of — along the plate. 

Ur\ 
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TABLE I 
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TABLE II 


VALUES OF 



FOR n = 7 


M 

( 

u l/ u i) 2 for 1 

= 7 and R & = 

— 



6 x 10 h 

1 X 10 5 

5 x 10 5 

1 x 10 6 

1.5 x 10 6 




■^w “ -^adw _ Q 






To 






0.2401 

0.1992 

0.1343 

0.1122 

0.1019 



.2538 

.2110 

.1427 

.1191 

.1086 



.2878 

.£412 

.1681 

.1408 

.1264 

3 

.4013 

.3310 

•2793 

.1955 

.1650 

.1492 

4 

.4492 

.3768 

.3204 

.2241 

.1903 

.1740 

5 

.5020 

.4216 

.3622 

.2528 

.2174 

.1992 




T w - T adv Q ^ 






To 




0 

0.3254 

0.2642 

0.2189 

0.1502 

0.1267 

0.1155 

1 

• 3384 

. .2757 

.2309 

.1576 

.1332 

.1222 

2 

.3721 

.3051 

.2567 

.1764 

.1498 

.1362 

3 

.4159 

•3443 

.2916 

.2034 

.1722 

.1570 

4 

.4629 

.3869 

.3297 

.2314 

.1971 

.1799 

5 

• 5084 

.4290 

.3683 

.2611 

.2228 

.2048 




Tv “ T a dv 

= -0.5 






To 




0 

0.2642 

0.2065 

0.1710 

0.1139 

0.0933 

0.0839 

1 

.2845 

.2266 

.1866 

.1236 

.1027 

.0931' 

2 

•3317 

.2683 

.2233 

.1502 

.1262 

.1144 

3 /'• 

.3851 

.3160 

.2669 

.1828 

.1548 

.1407 

4 - 

.4421 

.3652 

• 3114 

.2151 

.1834 

.1677 

5 . 

.4965 

.4147 

• 3544 

.2490 

.2125 

.1945 




Tv “ Tadv 





' 

To 

- = 1.0 



0 

0.3469 

0.2841 

0.2386 

0.1636 

0.1399 

0.1268 

1 

.3588 

.2940 

.2472 

.1697 

.1441 

.1317 

2 

.3831 

• 3193 

.2704 

.1875 

.1576 

.1450 

3 

.4280 

• 3556 

.3024 

.2107 

.1798 

.1624 

4 

.4701 

.3945 

.3401 

.2372 

.2025 

.1856 

5 

.5162 

.4369 

.3838 

.2663 

.2280 

.2084 




T w - .T a ^ w 

T 0 

= -1.0 



0 

0.2125 

0.1608 

0.1267 

0.0759 

0.06076 

0.05382 

l 

.2440 

.1892 

.1525 

.0961 

.07919 

.07076 

2 

• 3058 

.2460 

.2016 

.1329 

.1113 

.1002 

3 

• 3703 

.3020 

.2515 

.1712 

.1439 

.1310 

4 

.4305 

• 3558 

• 3020 

.2078 

.1766 

.1612 

5 

• 4855 

• 4077 

.3455 

.2426 

.2075 

.1892 
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TABLE III 

MEASURED QUANTITIES USED IN COMPARISON 


Point 

(See fig. l) 

M 

Tw “ T a cLvj 
op 

T v “ ^adv 
To 

% 

Ng 

1 

1.59 

. 67 

0.145 

0.385 x 10 6 

270 

2 

l.6i 

58 

.126 

.6 97 

422 

3 

2.15 

-153 

-.320 

• 585 

388 

4 

2.19 

-171 

-.365 

1-39 

799 

5 

2.52 

-376 

-.776 

• 7 

450 

6 

2.58 

-394 

-.837 

1-973 

1130 

T 

2.60 

49 

• 135 

.208 

180 

8 

3.12 

-90 

-.258 

.457 

331 

9 

3.60 

-509 

-1.312 

• 921 

561 

10 

3.69 

-745 

-1.842 

1.07 

634 
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Figure *1.- Comparison of experimental results (refs. 3 an< l M an <l 

present method. 
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